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The p-n junction has provided the basis for the semiconductor-device industry. Investigations of 
p-n junctions based on Mott insulators is still in its infancy. Layered Mott insulators, such as the 
cuprates or other transition metal-oxides, present a special challenge since strong in-plane correla- 
tions are important. Here we model the planes carefully using plaquette Cellular Dynamical Mean 
Field Theory with an exact diagonalization solver. The energy associated with inter-plane hopping 
is neglected compared with the long-range Coulomb interaction that we treat in the Hartree-Fock 
approximation. Within this new approach, "Dynamical Layer Theory" , the charge redistribution is 
obtained at the final step from minimization of a function of the layer fillings. A simple analytical 
description of the solution, in the spirit of Thomas-Fermi theory, reproduces quite accurately the 
numerical results. Various interesting charge reconstructions can be obtained by varying the Fermi 
energy differences between both sides of the junction. One can even obtain quasi-two dimensional 
charge carriers at the interface, in the middle of a Mott insulating layer. The density of states as a 
function of position does not follow the simple band bending picture of semiconductors. 

PACS numbers: 74.78.Fk, 71.30.+h, 71.10.Fd 



I. INTRODUCTION 

Strong electronic correlations in transition-metal 
oxides manifest themselves in many spectacular ways, in- 
cluding high-temperature superconductors, Mott insula- 
tors, exotic magnetic phases and colossal magnetoresis- 
tance materials. 1 . Although these phases are present in 
the bulk, we expect to obtain even more complex and in- 
teresting states of matter if we harvest the properties of 
correlated electron with the help of heterostructures 2-6 . 
A decade ago, Ohtomo et al. 7 demonstrated the technical 
capability to create atomically sharp interfaces between 
transition metal oxides, launching a revolution in the field 
of oxide heterostructures. 

The field is presently dominated by the well known 
LaAlOa/SrTiOs interface, which gained much popularity 
because it harbours a two-dimensional electron gas 8 . In 
the bulk, both oxides are insulators and non-magnetic 
materials. However, when we consider the interface, it 
can become superconducting 9 , ferromagnetic 10 and both 
phases can even coexist 11 . These effects can all be ex- 
plained by charge transfer. It is thus possible with an 
external applied electric field to explore the phase dia- 
gram 12-14 and create exotic interface effects 15 . 

While this much celebrated interface is a central theme, 
it is not the only interface currently under study. Charge 
transfer at the interface between non-superconducting 
cuprates was also observed, giving rise to interfacial su- 
perconductivity 16,17 . 

Different effects caused by electron transfer at the 
interface between a hole-doped Mott insulator and an 
electron-doped Mott insulator have already been obser- 
ved in cuprate heterostructures 16 ' 18 and the charge redis- 
tribution has been the subject of a few models 19,20 . There 
have been device proposals based on charge redistribu- 



tion in heterostructures involving Mott insulators 21 . 
Also, Mott insulating p-n junctions have been sugges- 
ted to potentially yield a photovoltaic effect with high 
energy-efficiency 24 . There is some evidence that two over- 
doped cuprates, one electron-doped, one hole-doped, ar- 
ranged in a p-n junction geometry, will yield a high- 
temperature superconducting effect due to the artificial 
doping provided by charge redistribution 25 . 

Inspired by the experiments on cuprates, here we focus 
on the electronic charge density distribution and local 
density of states at the interface between a layered hole- 
doped Mott insulator and a layered electron-doped Mott 
insulator. 

Theoretical understanding of such interface problems 
must first address the question of bulk strongly correlated 
materials that can exhibit Mott insulating behavior 26 . 
This requires advanced methods, such as dynamical mean 
field theory (DMFT) 27 . In low dimension, when the self- 
energy can acquire strong momentum dependence - as 
happens for example with d-wave high-temperature su- 
perconductors - cluster extensions of DMFT such as Cel- 
lular Dynamical Mean-Field Theory (CDMFT) or Dy- 
namic Cluster Approximation (DCA) are necessary. 28-33 
Potthoff and Nolting 34 initiated research on surfaces and 
interfaces with DMFT by introducing " inhomogeneous" 
layered DMFT, that was subsequently applied to various 
interfaces. 35-37 

The case of a Mott/conventional insulator junction 
was studied in Hartree-Fock theory (HFT) 38,39 and 
(DMFT) 40,41 by Okamoto and Millis. Similarly, the 
charge profile of the manganite/cuprate junction 42 and 
the polar/non-polar Mott insulator junction was des- 
cribed with both HFT and DMFT. Magnetic ordering at 
the interface and the corresponding rich phase diagram 
can also be explored with simple HFT 38,39,42,44 . While 
both HFT and DMFT can be used to obtain an estimate 
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of the charge redistribution, also known as electronic re- 
construction, HFT leads to a simple band bending pic- 
ture while the local density of states as a function of the 
position obtained from DMFT is non-trivial 36,37 ' 40 ' 45 . 

Here we provide theoretical tools beyond DMFT that 
are appropriate for layered materials. We take into ac- 
count both the non-trivial non-local correlations within 
the individual layers due to strong short-range repulsion 
using Hubbard model, and the charge redistribution bet- 
ween layers that is dominated by long-range Coulomb re- 
pulsion. In-plane quantum fluctuations are included but 
correlations between the planes are treated at the mean- 
field (Hartree-Fock) level, neglecting the energy associa- 
ted to hopping between planes. We call this approach 
Dynamical Layer Theory (DLT). We focus here on the 
normal state but it will be clear that our approach is 
easily generalizable to broken-symmetry states. 

The model and method can be found in Sec. II. The 
charge density profiles obtained both numerically and 
analytically, along with density of states profiles can be 
found in Sec. III. We focus the discussion in Sec. IV on the 
domain of validity of the method and on the difference 
between Mott p-n junctions and ordinary semiconductor 
p-n junctions, concluding with a summary and sugges- 
tions for simple extensions of the method. 



II. MODEL AND METHOD 
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Figure 1: (a) Schematic view of the interface. Each plane 
corresponds to one of the CUO2 planes of a cuprate material. 
The two colors represent 2 different cuprates. (b) Schematic 
representation of the relevant energies for both materials. On 
the left is the electron doped PCCO and on the right the 
hole doped LSCO. Values are not to scale. Due to the Mott 
nature of these material, there is a gap inside each of the 
energy bands, but this is omitted here. It is addressed later. 
Both Fermi energies e/,i and e/,2 are negative since we take 
the vacuum as the zero energy state. Ae/ is the difference 
between the Fermi energies. 



Although we will consider more general cases, we 
are motivated by p-n junctions made of a thin film of 
Pi'2-(5 e Ce ( 5 i ,Cu04 (PCCO) deposited over a thin film of 
La 2 _5 h Sr,5 fc Cu04 (LSCO). 25 Fig. 1 presents a schematic 
view of this p-n junction with the corresponding defini- 
tions of the Fermi energies and of their difference for the 
bulk systems. 

The Hamiltonian for this p-n junction is given by 

H = H ll +H ± + H c (I) 

where H± stands for the hopping between planes, He 
for the Coulomb interaction between planes and Hi is 
the Hubbard Hamiltonian of the planes 

"ii = £*f 

k 

with cfy the annihilation (creation) operator at site i 
in plane k for spin a and rv% a the corresponding number 
operator. The in-plane hopping matrix t^, Hubbard in- 
teraction U k and site energy Eq can depend on the mate- 
rial. Here we assume that ty for second neighbor hopping 
is t' = -O.m 46 , t" = 0.08t for third neighbor hopping 
and U = 8t independent of k while the site energy £q 
is chosen such that we obtain the required bulk filling 



in plane k. Note that in reality, U is somewhat smal- 
ler for the electron-doped cuprates when compared with 
hole-doped ones, 47,48 but taking the same value suffices 
to illustrate the physics. 

Our main approximation is to neglect H± in the Hamil- 
tonian Eq.(l) compared with H\\ and He- The large ani- 
sotropy in the hopping amplitudes t± <C t of the cuprates 
motivates this approximation. The physics of charge re- 
distribution will be correctly taken into account when the 
inequality (H±) (H c ) is satisfied by the final solution. 

The Coulomb interaction coming from the deviation 
from the bulk charge for the Vth plane is given by 

H c = ^ £ ^(2||) n m(-g||)Vi, m (cj<||). (3) 

0,771, (J II 

where ni(q^) is the in-plane Fourier transform of the 
charge-density operator for plane I and V^ m (q\\) the Cou- 
lomb potential. Self-interaction should not be included 
and it is understood that the uniform background must 
be subtracted from the charge-density operator for the 
qii = contribution to the sum. When charge redistri- 
bution occurs on large length scales, the Hartree-Fock 
approximation is justified for this piece of the Hamilto- 
nian. The Fock term being short range, we assume that 
it is taken into account by the Hubbard in-plane Hamil- 
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tonian. We are then left with 

HC ~ 2 ^ ( n i(<?||))«m(-9||)^,m(g||) 
l,m,q[\ 

~ \ Y ( n l(<l\\))( n rn{-q\\Wl, m (q\\) 



(4) 



where brackets denote thermal and quantum mechanical 
averages. 

While we could treat cases where the charge distribu- 
tion breaks lattice translational symmetry in the planes, 
here we assume uniform solutions. In other words only 
the qti = contribution to (ni(q\\)) survives. Because 
there are no terms left that involve quantum mechanical 
operators that couple different planes, the eigenstates are 
direct products of single-plane states and the grand po- 
tential obtained from (H), the ground-state expectation 
value of H, is 



(H)-HN = Y, 



E (rife) -I Ant - (m - e f,k)) n k 

(5)" 

where now stands for the expectation value of the 
qu = number operator in plane k that we rewrite expli- 
citly as Ant = rlfc— no,fc to emphasize the deviation, men- 
tioned above, with respect to the nominal charge hq^ in 
the plane. The quantity E(nk) is defined by (H^) and the 
chemical potential fi is a Lagrange multiplier that gua- 
ranties charge neutrality. (Note that if the Fermi energy 
e/,fe (see Fig.(l)) is independent of the plane k, then the 
chemical potential is equal to the Fermi energy and all 
the planes have their nominal bulk filling.) The electro- 
static potential energy ifk is given by the parallel-plate 
capacitor formula 



<ySfc (rti <k ) = - 2^2^ 

m=l 1 = 1 



k—X m / \ 

(ec,m + ec,m+i) 



An, (6) 



where we have set the zero of potential on the left-hand 
side of the junction and defined ec,k — e dk/-^-k e ki with 
e the fundamental charge, Ak the area of Cu02 in a unit 
cell, dk the distance between CuC>2 planes and e k the 
dielectric constant 49 . 

Our problem now reduces to minimizing the functional 
of the classical variables representing average occupation 
numbers rife. The energy E(rTk) is in general a non-linear 
function of that can be calculated exactly, or approxi- 
mately. It contains all the information about the in-plane 
quantum fluctuations. This is the approach we call Dy- 
namical Layer Theory (DLT). 

In two dimensions, the self-energy near half-filling can 
acquire non-trivial momentum dependence, even in the 
normal state. This is why the approximate method that 
we use to obtain E(rik) in the layers is CDMFT on a 



2x2 plaquette. Single-site DMFT, used in many pre- 
vious studies, would be inappropriate here. We use an 
exact diagonalization solver 32 ' 33 ' 50-52 , then we perform 
the minimization of the grand potential Eq. (5) to ob- 
tain the charge redistribution using both numerical and 
analytical methods. From this approach, we describe an 
emerging phenomenon near the interface : a Mott deple- 
tion plateau. One can obtain many other properties, in- 
cluding the local density of states profile associated with 
the charge redistribution and thus the analog of "band 
bending" of the device interface. Many similarities and 
differences with the classical semiconductor p-n junction 
are highlighted in the discussion. 



III. RESULTS 

In the first subsection below, we show the n(eo) re- 
lation and corresponding density of states for a single 
plane. Numerical and analytical results follow respecti- 
vely in the two subsequent subsections. 



A. Single-plane results 

Fig. 2 presents the results for the filling n as a func- 
tion of £o that plays the role of the single-plane chemical 
potential u = dE/dn in the case of a single plane. The re- 
sults are plotted in Fig. 2. In the lower panel, the filling 
n vs Eo clearly exhibits the incompressible Mott phase 
where n = 1 as long as eo is in the gap. The region of 
zero compressibility, called Mott plateau is characteristic 
of a Mott insulator. With the Green function obtained 
from the CDMFT calculation we can also compute the 
density of state (DOS) as a function of the chemical po- 
tential. It is shown in the top panel of Fig. 2. The DOS 
clearly exhibits the Mott gap present in the "band dia- 
gram" for every possible chemical potential. We see that 
the compressibility (nr 2 dn/d[i) is finite only when the 
Fermi level (w = 0) crosses a finite DOS. There is an im- 
portant feature in the DOS around half filling (n = 1). 
We note that unlike a traditional band insulator, the 
DOS is not constant when the chemical potential varies. 
For example, for electron doping some spectral weight of 
the lower Hubbard band transfers to the upper Hubbard 
band at the Mott insulator transition (eq ~ 2.5i, 5.5t). 



B. Numerical results for the p-n junction 

A simple numerical minimization algorithm for this 
Hamiltonian gives the filling (rife) for every plane k in 
the ground state, and thus the charge redistribution of 
the junction. Results for an an electron-doped system on 
the left side and a hole-doped one on the right side is 
shown in Fig. 3. The resulting profiles exhibit some si- 
milarities with the classical semiconductor p-n junction. 
Nominal doping (5 ,i = \n ,i - 1|), charging energy (e c ) 



Figure 2: (a) Color coded DOS as a function of single-plane 
chemical potential and energy. The Fermi level in the DOS 
is at to = (dashed line), (b) Corresponding filling (n) as a 
function of single-plane chemical potential (fi = £o here). The 
results are for parameters appropriate for LSCO/PCCO, in 
other words next nearest- neighbour hopping t' — — 0.17t 46 , 
second nearest-neighbour hopping t" = 0.08t and U = 8t. 
Energy units correspond to t = 1. We also take h = 1. 



and Fermi level difference (Ae/ = e/,2 — e /,i) have been 
chosen to reproduce profiles that highlight some simila- 
rities and differences between this Mott p-n junction and 
semiconductor junctions. The Fermi level is lower for the 
hole- than for the electron-doped Mott insulator so that 
charge is transferred from left to right. The Fermi level 
difference increases in absolute value from panels (a) to 
(c). 

As in semiconductors, there is charge exchange and 
some regions can become depleted. Unlike semiconduc- 
tors, where region b and d of Fig. 3(c) would extend all 
the way to the interface and thus be adjacent and form 
the well known depletion layer, here there is an additional 
region c in the middle. This region behaves like a two di- 
mensional electron gas (2DEG) in the middle of an Mott 
insulating depletion layer formed by b and d. Outside the 
layer formed by 6, c and d, there is Thomas-Fermi scree- 
ning in region a and e. Semiconductors could in principle 
be forced into this situation when there exists an appro- 
priate Fermi energy difference, but this does not happen 
in practice, as far as we know. 

Let us focus on the evolution of the charge redistribu- 
tion profile caused by the Fermi level difference (Ae/) in 
Fig. 3. The increase in charge transfer, caused by an in- 
creased Fermi level difference, translates the reconstruc- 
ted charge profile, revealing more of the 'rigid' curves on 
each side without affecting their respective shapes, an 
effect we will discuss in the following subsection. 

A more striking difference with semiconductors arises 
when one uses the results for the DOS as a function of 
filling (Top of Fig. 2) , and the filling as a function of the 
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Figure 3: Examples of charge redistribution : The charge 
profile Ant is plotted as a function of layer number k. The 
parameters are U = 8t, es = 0.05t, no,i = 1.20, no, 2 = 0.88 
and Ae/ = — 2t, —5.5t and — 8t for (a), (b) and (c) respec- 
tively. These parameters have been chosen to shed light on 
qualitatively different features of the electronic charge redis- 
tribution. Blue is associated to the e-doped Mott insulator 
and red is associated to p-doped Mott insulator. The dots are 
obtained from the numerical calculation and the shading from 
the analytically calculated profile with the same parameters 
within the continuum approximation. The lighter shading in- 
dicate the position of the Mott depletion plateaux. 



position to compute the local DOS as a function of po- 
sition. This analog of the semiconductor heterojunction 
band diagram is illustrated in Fig. 4 for another set of 
parameters. For this example, nominal doping on each 
side is chosen in order to obtain only one Mott deple- 
tion plateau on the left part of the junction. Some band 
bending is observed in the profile, but there is also a 
band modulation around the Mott transition region (on 
the edge of the Mott depletion plateau). This is caused 
by the spectral weight transfer associated to the Mott 
transition. 

Note that there is no need for a doped Mott insulator 
on the right-hand side of Fig. 4 to obtain a Mott de- 
pletion plateau on the left. The right-hand part of the 
structure could easily be a semiconductor or an insulator 
with appropriate relative Fermi level. As long as there is 
sufficient charge transfer, there will be this plateau. In 
other words, the resulting charge profile on one side is 
not only rigid, but also its shape is independent of the 
material on the other side of the junction that attracts 
(or repels) the charges. 
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Figure 4: (a) Charge redistribution as a function of plane 
number (the z-axis) of the structure (V = 8t, Ae/ = —It, 
ec = O.lt). (b) Color coded DOS in arbitrary units as a func- 
tion of uu and plane number. The Fermi level in the DOS is at 
ui = (dashed line). We observe an insulating region several 
planes wide, where n = 1. 



C. Analytical results for the p-n junction 

For a large number of C11O2 planes, a numerical ap- 
proach is needed to obtain the filling from the mini- 
misation of Eq. (5). Nevertheless, the profile can be 
estimated analytically by taking the continuum limit 
— ^ 2 f dz. In this limit, the minimization procedure 
consists in finding the charge distribution n (z) where the 
functional derivative of the ground state energy Eq. (5) 
with respect to the charge distribution n (z) vanishes. 
This leads to a Thomas-Fermi equation 



fj, 



e f . z + -Aiptot 



¥>0) 



(7) 



where we defined A/i (z) = SE (n(z)) /Sn(z) and Aip to t as 
the total electrostatic potential-energy difference between 
the right-hand side and the left-hand side of the junction. 
This term arises from the functional derivative of the 
electrostatic potential energy 



S(Jdz'ip(z')An(z r )) 
5n{z) 



2cp(z) - Acp tot . (8) 



In the limit where both sides of the junction are semi- 
infinite, we have 



Aiptot = e/,2 - e/,i 

A* = (e/,2 + e/,i) /2. 



(9) 
(10) 



The second derivative of the Thomas-Fermi equation 
Eq. (7) gives us an equivalent of the Poisson equation 



d 2 d 2 A[i{z) jdz 2 — —ecAn(A/j,(z)). The main source 
of simplification comes from the fact that the relation 
between An and A/i(z), which can be obtained from 
to the lower panel of Fig. 2, can be approximated by 
a piecewise continuous function that is either linear or 
constant in the vicinity of half-filling. 

The differential equation can then be solved separately 
in different segments with boundary conditions between 
the segments given by the continuity equation for the 
potential and the electric displacement. The analytical 
solution for the case illustrated in Fig. 3(c), for example, 
appears as the shaded area. In regions c, d and e on the 
right-hand side, the charge profile is a piecewise conti- 
nuous function given by 



An c (z) = 5 ,2 
An d (z) = S ,2 
An e (z) = <5 ,2 



~(z-*2) 



~(*-*3 
C X TF 



Az 



sinh 



(Z-Z2) 



X 



TF 



(ii) 



In the depletion region d, the electron density difference 
Arid is equal to the nominal hole doping of the material 
6h- The term g (z) — exp [{z — {2z/± — z^)) /Xtf] is ne- 
gligible for structures that are thick enough (24 — > 00). 
2^4 — Z3 is the mirror position of Z3 with respect to the 
right end of the structure. The difference between Z3 and 
Z2, the length of the Mott depletion plateau, is given by 



Az = \tfWi + 2\E g /5^ 2 \ - 1) 



(12) 



where Eq is the Mott energy gap, Sfii = 5o y idn/dn the 
chemical potential difference between the nominal filling 
and half- filling, while Xtf = d/ ' \fec dn/dfi is the Tho- 
mas Fermi length-scale for the charge profile, all in the 
approximation dn/ d[i\ n _ Ll constant outside the Mott pla- 
teau. The result for the size of the Mott plateau Az, 
comes physically from the fact that there is an electric 
field that changes linearly with distance inside the pla- 
teau leading to an electrostatic potential difference that 
grows quadratically with distance. When the electrosta- 
tic potential energy difference becomes equal to the Mott 
energy gap Eg, the plateau ends. 

The function g (z) in Eq. (11) allows the electric field 
to vanish at the surface of the junction. Here, we want 
to stress the idea that, except for g (z), the shape of the 
piecewise continuous function describing the electronic 
density profile is rigid and does not deform with charge 
transfer : A larger Ae/ only reveals more of the curve 
(Eqs. (11))- The only remaining unknowns for the elec- 
tronic density profile is the position z 2 (and hence z% 
that just follows from Az) and z\ relative to the interface. 
These can be found from the continuity of the electric dis- 
placement at the interface and from A/i(0~) — A/i(0 + ) = 
Ae/ = e/.2 — e^i that follows from the Thomas- Fermi 
Eq. (7). Note that the energy scale ec is very difficult to 
obtain from first principles. It will vary considerably de- 
pending on material. Its value changes considerably the 
length-scale of the charge profile. 
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IV. DISCUSSION 

The Mott depletion plateau obtained here is consistent 
with the result obtained by Oka and Nagaosa 53 even 
though these authors used a completely different tech- 
nique, namely the one-dimensional density-matrix renor- 
malization group. In other words, our method, which ne- 
glects hopping energy between planes, predicts the same 
incompressible state as a full one dimensional treatment. 
This plateau is also present in some HFT 42 and DMFT 34 
studies. 

We can verify the validity of our approach (DLT) in the 
vicinity of the Mott plateau by comparing the hopping 
energy t± with the potential energy difference between 
two adjacent planes in the plateau ~ d dAfi(z)/dz. We 
can estimate the latter from the smallest value of the 
derivative, namely at the beginning of the Mott plateau. 
Since the profile just before the plateau is exponential, 
we find that 

t± < 8fii(di/\ T F,i) = S .iec,i(XTF.i/di) (13) 

must be satisfied for both sides i of the p-n junction. 
Since Xtf,i ^ di, ec,% > t and t > ij. in cuprates, 
this condition will hold even for low dopings Sq^. Ho- 
wever, we see that for nominal dopings 5gi that are 
too small, this condition will not hold. Indeed, we see 
from the expression for the size of the Mott plateau, 
Eq.(12), that small nominal dopings 5o,i imply a very 
large Mott plateau, consistent with small charge transfer 
and hence small long-range Coulomb interaction. Note 
that in the examples used in this work, the relative die- 
lectric constant is much larger than that which is usually 
quoted for the cuprates (~ 10 to 20). This leads to much 
larger values of the Thomas- Fermi screening length than 
the value Xtf,z ~ <k expected in the cuprates. 17 

There is an important parallel that can be drawn bet- 
ween this Mott p-n junction and the classical semicon- 
ductor p-n junction. The charge profile at the interface 
in middle panel of Fig. 3 corresponds to the depletion ap- 
proximation in semiconductors. There are however three 
major differences. 

1. The Mott gap can be crossed at the interface and a 
two-dimensional gas of charge carrier can appear in the 
middle of a depletion layer like region c in Fig. 3. Using 
the equation for the size of the depletion layer Eq. (12) 
but for semiconductors, the doping energy <5/x 2 is so small 
compared with the gap energy Eq that the length of the 
depletion plateau becomes very large. 

2. In a Mott p-n junction, the length-scale for the 
charge redistribution is much smaller because the char- 
ging energy ec is generally much higher. 

3. The most important difference is that in addition to 
band bending, there is also band modulation around the 



Mott depletion plateau region in contrast with rigid band 
bending in semiconductors. This might lead to a broader 
spectrum of light absorption and emission. 

We emphasize that the whole analysis has treated ec 
and Aef as free parameters. For PCCO and LSCO, the 
parameter ec should be around ten times larger than 
what we have used. The resulting profile would be ten 
times shorter that what we showed here. We do not know 
precisely the value of Ae t , but it should be of the order 
of the kinetic energy ~ t. We need a Ae/ larger than the 
Mott gap Ec to obtain a two-dimensional gas of charge 
carrier. This might be achieved through the application 
of a large potential difference. Overall, our results are not 
limited to cuprates but apply to all layered doped Mott 
insulators. 



V. CONCLUSION 

In summary, we developed a very simple and intuitive 
method, Dynamical Layer Theory, to account for Hub- 
bard model correlations and calculate the charge profile, 
local DOS profile and other observables of layered cor- 
related systems. The main approximation is to neglect 
interlayer hopping compared with the charging energy 
treated in the Hartree-Fock approximation. Once the so- 
lution is found as a function of chemical potential in the 
various planes, we are left with a minimization problem 
to find the filling of each plane. Here we found the solu- 
tion of the Hubbard model within the planes using pla- 
quette CDMFT with an exact diagonalization solver. We 
neglected broken symmetries and finite temperature ef- 
fects, but they can all easily be included within our ap- 
proach. This opens, for example, the possibility to treat 
d-wave superconductivity accurately in such structures. 

Both numerical an analytic calculations predict an in- 
compressible Mott depletion layer analog to the semicon- 
ductor depletion layer, and the possibility to obtain quasi 
two-dimensional charge carriers of opposite signs on ei- 
ther side of the interface. The calculated density of states 
near the Mott depletion plateau deviates from expecta- 
tions from simple band bending ideas and should have 
interesting consequences for the absorption and emission 
spectrum of this device. 

We are indebted to Patrick Fournier, Maxime Dion 
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was partially supported by a Vanier Scholarship from 
NSERC (M.C.) and by the Tier I Canada Research Chair 
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